Interplay between exchange interactions and charging effects in metalhc grains. 
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We study the effect of the exchange interactions in small superconducting grains near the Coulomb 
blockade regime. We extend the standard description of the grain in terms of a single collective 
variable, the charge and its conjugated phase, to include the spin degree of freedom. The suppression 
of spin fluctuations enhances the tendency towards Coulomb blockade. The effective charging energy 
and conductance are calculated numericaUy in the regime of large grain-lead coupling. 
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Introduction. Coulomb blockade in metallic grains is a 
well studied phenomenon 1,2]. The transport through a 
grain in the Coulomb blockade regime can be studied us- 
ingrate equations when the coupling to the leads is weak 
[l|, ■ The renormalization of the charging energy when 
the coupUng to the leads is large is also well understood 
0, 0, la ■ The strongly coupled regime is best stud- 
ied by introducing a single collective degree of freedom, 
the phase, conjugated to the total number of electrons 
in the grain Q/e ^ St M- The use of this variable is 
justified when the separation between the electronic lev- 
els within the grain can be neglected, or, alternatively, 
when the conductance of the grain is large. In this limit, 
the interaction effects within the grain can be described 
by a simple hamiltonian T^, expressed in terms of the 
total charge {Ec{Q — Qo)^), the total spin (JsS^) and 
the individual electronic degrees of freedom. In the pres- 
ence of attractive interactions in the grain a pairing term 
i^BCsT^T) should also be included ^ that wiU drive 
the system towards superconducting state with energy 
gap A. In such case both the pairing Xbcs and the ex- 
change Js will grow under renormalization group (inte- 
gration of energies down from Thouless energy) to a scale 
much larger than the bare one, which is initially of the 
order of the level spacing. Moreover due to the attrac- 
tive character of interaction the spin susceptibility due 
to exchange will be positive, so that spin fluctuations get 
suppressed much in the same way as charge fluctuations 
do due to the charging energy. The essential distinction 
between the two will stem from the topological differ- 
ences between spin group SU(2) and charge group U(l) 
in which their conjugate phases exist. 

In the following, we will generalize the usual descrip- 
tion of a small grain in the Coulomb blockade regime in 
terms of phase dynamics 0- We include also the dy- 
namics of the total spin of the grain in SU(2), on the 
same footing as the total charge. We assume that the 
grain has a negligible level spacing and a finite positive 
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renormalized susceptibility according to the arguments 
above. 

The coupling between the grain and an external (nor- 
mal) electrode leads to single electron tunnelling. We 
describe it as usually done, in terms of a long range in- 
teraction, in time, between the phase and spin variables 
of the grain In most of the calculations, we neglect 
the long time cutoff in the kernel describing these inter- 
actions imposed by the superconducting gap. This ap- 
proximation is reasonable when the bare charging energy 
is larger than the gap, and its renormalization arises from 
virtual tunnelling processes of high energy. The effects 
arising from the suppression of the subgap conductance 
due to the superconducting state can be incorporated 
into our calculations in a straightforward way. Finally, 
our method can be used to study superconducting junc- 
tions with subgap leakage currents [l^ , although wc will 
not consider this case in detail. 

The effects of a constant exchange term on the trans- 
port properties of a quantum dot has already been stud- 
ied in the limit where the coupling to the leads is weak, 
using rate equations [13 . ITEl fl^ . The present formal- 
ism goes one step beyond this by summing all processes 
up to cotunneling level. However we do not consider here 
the changes in the grain susceptibility induced by the 
spin-orbit coupling which has been considered by other 
authors [13 . 

The model. As mentioned above, we will focus on the 
case in which the superconducting gap is smaller or of 
similar order as the renormalized charging energy, E^, 
and analyze how this renormalization of the bare charg- 
ing energy, Ec^ depends on the exchange J. These con- 
straints are satisfied, for instance, by Al superconducting 
grains with radii below 100 nm. 

The hamiltonian that we study is: 7Y Hgrain+^ioad+ 
Hhop, where 

k,s 

Hhop = -tJ2cldk,s + h-C. (1) 
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and N and S are the total number of electrons and 
the total spin of the grain, TV = J^i s ^1 s'^i s ^^'^ ^ = 

5 Si s s' ^1 s'^s s''^i s'' where a denotes the Pauli matri- 
ces. The grain-dot conductance, in dimensionless units, 
can be approximated by a w i^Pgrain(eF)picad(eF), where 
Pgrain/ioad(eF) is the density of states at the Fermi level 
of the grain and the lead. As mentioned above, we ne- 
glect the energy dependence of the density of states of 
the grain. The only interactions included in eq.(^ are 
through the total spin and charge of the grain. 

Path integral formulation. We can integrate out the 
fermionic degrees of freedom and obtain a description 
in terms of collective variables only by using the path 



integral formalism. The action is 5 = ■S'^rain + "^i 
Shop + Sint, where 



lead 



9" 



Si 



lead 



di,s [dr + Ci - //grain) d^,s 
I dr'S^ Ck.s {dr + Cfc - Mload) Ck,s 

-t / drVd 



^hop = -t I dr y di^s Ck,s + H.c. 

5int = EciN ~ N,^tf + JsiS - S,^tf (2) 

We have included an offset electron number iVoxt = 
V^ato /eCg induced by a gate voltage V^atc and an offset 
spin 5'oxt = -ffcxt/2Js induced by an external magnetic 
field Hcxt, that couples to the total spin. 

We can now decouple the quartic interaction term 
Sint by means of a Hubbard-Stratonovich transformation 
e-Ec{N-N^^,f oc /pv^e-^'/4Bc-»v(JV-iv„0, and simi- 

larly for 5, which introduces a new scalar field V for the 
total charge and a vector field H for the total spin. We 
then have iSint — Sq + Si, with 



So 



Si 



l3 

dr 



iVNc; 



iH ■ Sc_ 



dT{VN + H-S 



(3) 



We now perform a time dependent canonical transforma- 
tion (a phase and spin rotation) on the electronic wave- 
functions, in order to cancel the term Si in eq.®. This 
U{1) X SU{2) transformation can be written as: 



dksir) 
U{t) 



Uss'It) dks'ij) 

e 



(4) 



The transformation is parametrized by the angles 4'{t) 
and ^(r), and by the three dimensional unitary vector 
n(T). The requirement that 5*1 in © is cancelled implies: 



[drU) =iV+-H-a 



(5) 



so that V = (j) and 

H = S_h + sin ^ n + (1 — cos f ) 



n X n 



(6) 



These identities provide an alternative and convenient 
parametrization of the auxiliary fields V and H, repre- 
sented now by the (f>, ^, n fields, which will be used in 
the following. Note that the U{1) gauge transformation 
needed to replace V by the phase </> leads to the standard 
description of charging effects in terms of phase fluctu- 
ations. It is interesting to note that eq.® implies that 
H is proportional to the angular momentum of a sphere, 
considered as a rigid body . The periodicity in imag- 
inary time of the arguments in the action implies that 
U{0) — U{(3). This constraint implies the usual quanti- 
zation of the charge in the grain, and also of the spin (see 
below), due to the discreteness of transport events. 

A more compact notation for the transformation in 
eq.Q can be given in terms of the following, r dependent, 
two- and four-dimensional unit vectors: 



u-r ~ (sin cos 0r) 



fir sin i^^, cos 



(7) 



We can use these vectors to write Sq as: 

where external gates and fields have been taken as zero 
for the moment. The transformation in eq.|0J modifies 
also the lead-grain coupling. Shop'- 

Shop ^ dr Ci^s Us,s' dk^s' + H.c. (9) 

Jo , , 

A final step to obtain the effective action for the ro- 
tor fields is to integrate out the fermionic fields to order 



- p-|(sLp>o+o(t") 



t^, using (e '^''°p)o = e ^i^ha-p 
following dissipation term: 



One obtains the 



S'di: 



/5 

dr 



dr' K{t - t') Tr [c/^^C/^, + ulu^ 



where /^(t) = -[G'lead(T)Ggrain(-T)]/[picad(eF)Pgrain(eF)] 

(ttT)^/ sin^(7rTT), Giead and Ggrain being the lead and 
grain unperturbed Green's functions in imaginary time. 
Recall here the that the finite range that the supercon- 
ducting gap could bring in is assumed larger than the 
decay time of the phase correlators which is of order , 
so that the gapless K{t — t') of the normal state yields 
equivalent results. S'diss may be finally recast as 

•S'diss = dr dr' K{t — t') 1 — {dr-UT') {vt 'Vt') (10) 
Jo Jo 

This term is sufficient to account for second order tun- 
nelling processes, and in particular it can describe cotun- 
nelling features. The derivation is valid when the conduc- 
tance between the grain and the electrode per channel is 
small, and it can be used even if the total value is large. 



3 



The method leading to ea. H10() can be easily general- 
ized to the case when the density of states in the grains 
or in the leads is spin dependent. The effective action 
will contain terms involving sin(^T- +^t') which break the 
symmetry between the four components of the vector Vt- 
These terms are analogous to the Josephson term which 
arises in the charge dynamics when the leads, or the grain 
are superconductors 7]. 

The final action is Scs = 5*0 + S'diss: written in terms 
of the dynamical variables u,- and Vr only. In the limit 
Js = 0, the field Vr can be taken as a constant, and the 
model reduces to the standard phase only model. 

Results. It is instructive to analyze first the decou- 
pled grain, described by 5*0 in eq.®, to see where this 
spherical rotor description of the total spin comes from. 
As mentioned above, 5*0 contains the usual phase term, 
which leads to the quantization of the charge, and a con- 
tribution which is equivalent to that of a rigid rotor, 
and which leads to the conservation of spin. The eigen- 
values associated to 5*0 can be written as En,s.s,,k — 
EcN'^ + JsS{S + l), where 0, 1, 2 • • • , 5 = 6, 1,^2 • • • , 
< S:, < S and -S < K < S. The degeneracy of a 
given state is (25* -I- 1)^ (l^. This degeneracy can be 
understood by noting that, in the limit studied here, the 
level spacing within the grain can be neglected. The grain 
energy is solely determined by the total charge and the 
total spin. Let us assume that, in the neutral dot, there 
are Nq spin 1/2 electrons which contribute to the total 
spin. The number of states of total spin S (each with 
degeneracy 25* -f- 1) is: 



No 



No 



(11) 



In the limit of many electrons Nq — > oo. 



lini 



N„/S- 



= (25 + 1)Cno, where C, 



No 



one obtains 

_ 2"0 + 3/2 . 



a constant independent of S. This means that the total 
degeneracy of a state composed of many 1/2 spins and 
given value of the total spin momentum (S"^) = S{S + 1) 
is CnoC^S +1)'^ , just as Cn^ rigid rotors with total angu- 
lar momentum S. The existence of this degeneracy leads 
to a prefactor in the free energy which is independent 
of the angular momentum. This multiplicity, like similar 
degeneracies in the case of ordinary Coulomb blockade, 
does not affect the effects associated to the spin gap dis- 
cussed in this paper. 

The following calculations including the full action S'off 
have been done by averaging over all paths in the unit cir- 
cle parameterized by u and the four-dimensional sphere 
which defines v, using an extension of the Monte Carlo 
code developed earlier for related problems 0, . The 
effective charging energy is calculated by summing over 
winding numbers of the phase. The conductance between 
the grain and the electrode has been approximated by 
the expression G{(3/2) 0, |23|, vahd at low tempera- 
tures, where G is the correlation function, in imaginary 
time, of the variable VtUt- The latter combination de- 
scribes the transfer of a full electron to the grain. We 
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FIG. 1: Charge-charge, spin-spin and electron-electron cur- 
rent correlations (see text) versus inverse temperature for 
a = 0.3 and different values of Js- 
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FIG. 2: Renormalized charging energy of the grain in the pres- 
ence of a finite spin gap Js, versus the dimensionless grain- 
lead coupling a. Note that the decay becomes less pronounced 
for growing spin gap. 



calculate, separately, the correlations Gu = {ut 'Ut') and 
Gy = {vr ■ Vr') which correspond to charge only and spin 
only currents. 

The current correlation functions are shown in Fig. . 
It is interesting to note that both {ur ■ Ut') and {vr ■ 
Vr') decay exponentially, while the composite correlation 
{ur ■ Ur' Vr ' Vr') decays as (r — t')~^, as required by Grif- 
fith's inequality [23|. The differences between the phase- 
phase, "rotation-rotation" and current-current correla- 
tions is reminiscent of the behavior of a Luttinger liquid. 
It implies that the electron current cannot be factorized 
into its spin and charge components. The exponential 
decay of the correlations associated with the collective 
charge and spin degrees of freedom can be understood as 
the effect of a charge and spin gap in the grain. It can 
be obtained by making a mean field decoupling of the 
variables, in a similar way to the calculation for charg- 
ing effects in coupled grains The (r — r')~^ decay of 
the current-current correlation describes the cotunnelling 
processes at low temperatures. 

The effective charging energies, as functions of a and 
Js, are shown in Figs. [2] and [3]. The effect of a finite Js 
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FIG. 3: Renormalized charging energy of the grain in the 
presence of a finite spin gap Js, versus the value of the spin 
gap Js- Note the saturation for large Js- 



on the renormalized charging energy is significant, even 
for small values of Js- We can estimate analytically this 
effect, by assumming that when Js — > the fluctuations 
in the variable Vr are small. The effect of these fluctua- 
tions on the variable Ut can be approximated by replacing 
a in ea. (|10|) by a (j-Dp). Assuming that the fluctuations 
of Vr are harmonic, we find: 



Ec 



(12) 

ujy2Js Ec ^ ' 



where A is a high energy cutoff, comparable to the elec- 
tronic bandwidth. Then, using the well known expression 
for the renormalized charging energy for large values of 

a H 0: exp ^-2Tr^a - |. This en- 

hancement of the effective charging energy by a spin gap 
is another manifestation of the non-separability of charge 
and spin. 



Conclusions - We have analyzed the influence of the 
exchange term in a small superconducting grain on the 
charging effects in the regime where the superconduct- 
ing gap is smaller or comparable to the charging energy. 
The suppression of the spin susceptibility reduces large 
fluctuations in the spin of the grain, and enhances the 
tendency towards Coulomb blockade. Our analysis inte- 
grates out the electronic degrees of freedom in the grain 
and in the external leads, and provides a simple descrip- 
tion in terms of the charge and spin degrees of freedom 
of the grain only. 

The effects of the exchange term have been analyzed 
for closed quantum dots , w hich are almost decoupled 
from the leads 0, lisl IT^ . Our scheme provides a 
generalization which is non-perturbative in the coupling 
strength in the sense that one can recover exponential 
effects in the coupling, such as the renormalization of 
the charging energy, which cannot be derived from the 
addition of sequential processes. 

A statistical approximation to the electron-electron in- 
teractions in a small dot predicts that the bare exchange 
term is negative and of the order of the separation be- 
tween electronic levels [13, El- Our analysis, on the other 
hand, is valid only when the exchange term is positive 
and larger than the level spacing. This regime corre- 
sponds to systems with an attractive electron-electron 
interaction near a superconducting transition, when the 
exchange J is significantly enhanced 24] . Spin fluctua- 
tions in a superconducting grain at low temperatures can 
therefore have a strong influence on charge fluctuations, 
restoring the system to a Coulomb blockade regime even 
when the coupling to the leads is strong. 
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